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ABSTRACT 

Direct measurements of dark matter distributions in galaxies are currently only possible through the 
use of gravitational lensing observations. Combinations of lens modelling and stellar velocity dispersion 
measurements provide the best constraints on dark matter distributions in individual galaxies, however 
they can be quite complex. In this paper, we use observations and simulations of gravitational 
microlensing to measure the smooth (dark) matter mass fraction at the position of lensed images 
in three lens galaxies: MG 0414+0534, SDSS J0924+0219 and Q2237+0305. The first two systems 
consist of early- type lens galaxies, and both display a flux ratio anomaly in their close image pair. 
Anomalies such as these suggest a high smooth matter percentage is likely, and indeed we prefer 
~ 50 per cent smooth matter in MG 0414+0534, and ~ 80 per cent in SDSS J0924+0219 at the 
projected locations of the lensed images. Q2237+0305 differs somewhat in that its lensed images lie 
in the central kiloparsec of the barred spiral lens galaxy, where we expect stars to dominate the mass 
distribution. In this system, we find a smooth matter percentage that is consistent with zero. 
Subject headings: dark matter — gravitational lensing — quasars: individual: MG 0414+0534 - 
quasars: individual: SDSS J0924+0219 — quasars: individual: Q2237+0305. 



1. INTRODUCTION 

Gravitational lensing measurements provide the only 
direct method to probe the non-luminous matter com- 
ponent of lensing systems. Weak lensing measurements 
of cluster masses and mass distributions are becoming 
routine, however the measurement of the dark matter 
components of individual galaxies is a new and exciting 
prospect. By analysing gravitational microlensing sig- 
nals in some strongly lensed quasars, we are able to infer 
the ratio of clumpy to continuously distributed matter 
along the line of sight to the background source. This 
gives us a probe of the dark matter content in a lensing 
galaxy at p ro jected radii of ~ 2-10kpc from the centre of 
the g alaxy (jSchechter fc Wambsganssl l 200llPoolev et all 

Recent years have seen considerable advancement 
in the mapping of dark and stellar mass in lens- 
ing galaxies (e.g. iKeeton. Kochanek. fc Falcol 119981 : 
iFerreras. Saha fc Williamsl 120051: IBarnabie et all 12009) . 
Usually, the total mass within the Einstein Radius is 
constrained by modelling the lensing galaxy to fit the 
observed lensed image positions. Photometry of the lens- 
ing galaxy, in combination with stellar population syn- 
thesis models, provides the distribution of stellar mass. 
Finally, observations of stellar velocity dispersions can 
be used to break the mass-sheet and mass-anisotropy 
degeneracies, and thus constr ain the overall mass den- 
sity profile (see fo r exam pl e iKoopmans fc Treul 120031: 
Treu fc Koopmansl 12004 IFerreras. Saha. fc Williamsl 



2005 1 : IFerreras. Saha. fc Buries! 120081: IBarnabie et al.l 
2009: lAuger et alj|2009l and references therein). 
These analyses can provide a very detailed picture of 
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the structure of lensing galaxies. However, they are quite 
complex, relying on often difficult observations and de- 
tailed modelling. A complementary method exists, in 
which observations of microlensing in the lensed quasar 
images are used to constrain the dark matter percentage 
along those lines of sight directly. 

Cosmological microlensing occurs when the light path 
to a lensed quasar image intersects a starfield in a fore- 
ground lensing galaxy. The lensing galaxy as a whole 
magnifies the lensed image; microlensing by individual 
stars induces variations about this macro-magnification. 
This is most readily detected in lightcurves of quasar im- 
ages, where relative motion between observer, lens and 
source causes uncorrelated fluctuations in brightness be- 
tween lensed images. This ef fect was first obse rved in the 
lensed quasar Q2237+0305 (jlrwin et alJll989f ). 

In some lensing systems we observe a close pair of 
quasar images. Basic lensing theory suggests that these 
two ima ges should have approx i mately equal magnifi- 
catio ns (|Chang fc Refsdall Il979l : iBlandford fc Naravanl 
1986). On the contrary, in eight out of ten known cases 
we find that the brightness of the image located at the 
saddle point in the time delay surface is suppressed rel- 
ative to the image at the minimum in the time delay 
surface (|Poolev et al.1 1200 7). Microlensing is one possi- 
ble explanation for such flux ratio anomalies. However, 
micro lensing by a purely stellar comp onent is not suffi- 
cient. iSchechter fc Wambsgansa (|2002f ) showed that this 
discrepancy could be accounted for by adding a signifi- 
cant smooth matter component to the lens at the image 
positions, since minimum and saddle point images are 
microlensed differently when a smooth matter compo- 
nent is added. 

We have previously developed a technique for using 
single-epoch multi- wavelength observations of anomalous 
lensed quasars to constrain the radius and radial pro- 
file o f the background quasar accret ion discs (|Bate et al.1 
I2008t iFlovd. Bate, fc Webster! 12009ft . In those analyses, 
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we marginalised over the smooth matter percentage in 
the lens as a nuisance parameter. Here, we turn the 
problem around and instead marginalise over the quasar 
parameters to obtain constraints on the smooth matter 
percentage in the lens at the image positions. 

Rough microlensing measurements of smooth matter 
percentages have been reported prev iously. Spectroscopy 
of SDSS J0924+0219 undertaken bv lKeeton etH] (l2006f) 
suggested a smooth matter percentage of 80 to 85 per 
cent in that lens at the location of the D and A images. 
Using X-ray monitoring of HE 1104+1805. lChartas et al.l 
(2009) reported a s mooth matter perce ntage of ~ 80 
per cent is favoured. iPoolev et al.l (|2009l ) measured the 
smooth matter percentage in PG 1115+080 to b e ~ 90 
per cent, using X-ray observations. iMorgan et al.l (2008) 
found a w e ak tre nd supporting this result. Most recently, 
iDai et afl ()2009h favoured a smooth matter fraction of 
~ 70 per cent using X-ray and optical monitoring of 
RXJ 1131-1231. Microlensing analyses consistently pre- 
dict a significant smooth matter percentage in the lensing 
galaxy at the position of anomalous images. 

In this paper, we present constraints on the dark 
matter percentages in three lensing galaxies: MG 
0414+0534, SDSS J0924+0219 and Q2237+0305. MG 
0414+0534 and SDSS J0924+0219 are both lensed by 
early-type galaxies, and consist of close image pairs dis- 
playing a flux ratio anomaly. MG 0414+0534 is moder- 
ately anomalous, whereas SDSS J0924+0219 is the most 
anomalous lensed quasar currently known. Q2237+0305 
differs from the previous sources in two key ways: it is 
lensed by a barred spiral galaxy, and it does not contain a 
close image pair. Nev ertheless, it is kno wn to be affected 
by microlensing (e.g. Ilrwin etanH989h . 

This paper is laid out as follows: in Section[5]we discuss 
the observational data on the three systems of interest. 
The simulation technique is briefly described in Section 
[21 We present our results and discussion in Section @] 
and conclude in Section [SJ Throughout this paper we 
use a cosmology with Ho = 70kms~ 1 Mpc _1 , Vt m = 0.3 
and Vl A = 0.7. 

2. OBSERVATIONAL DATA 
2.1. MG 04U+0534 

MG 0414+0534 was discovered bv lHewitt etail (|1992t ). 
It consists of a bac kground quasar at z s — 2.64 
(|Lawrence et al.lll995D a nd a foreground early-typ e lens- 
ing galaxy at zi = 0.96 (|Tonrv k Kochaneldll999h . Four 
images of the quasar are observed, with the close im- 
age pair (images A\ and A2) displaying a flux ratio 
anomaly. This anomaly is weak in both t he mid-infrared 
(A 2 /A 1 = 0.90+0.04 on 2005 October lO. lMinezaki et al.l 
2009) and the radio {An /A x = 0.90 + 0.02 on 1990 April 
2. iKatz k Hewitt! 119931 ). but somewhat stronger in the 
optical (A 2 /A 1 = 0.45 ± 0.06 on 1991 November 2-4, 
ISchechter k Moord[l993T) . 

In our analysis, we used three epochs of multi- 
wavelength MG 0414+0534 observations, presented in 
Table [1] The first two epochs were archival HST 
data, obtained from the CAS TLES Survey webpage0 
(jFalco. Lehar. k Shapiro! 119971 ) . The third epoch was 
obtained by us using the Magellan 6.5-metre Baade tele- 

3 http:/ /cfa-www. harvard.edu/castles/ 



Table 1 

Observed flux ratios in MG 0414+0534 



Filter 


A c (A) 


Fobs = A 2 /A 1 


Date 


H 


16500 


0.67 ±0.05 


2007 November 3 


.1 


12500 


0.60 ±0.2 


2007 November 3 


z' 


9134 


0.34 ±0.1 


2007 November 3 


V 


7625 


0.26 ±0.1 


2007 November 3 


r' 


6231 


0.21 ± 0.1 


2007 November 3 


F205W 


20650 


0.83 ± 0.03 


1997 August 14 


F110W 


11250 


0.64 ± 0.04 


1997 August 14 


F814W 


7940 


0.47 ±0.01 


1994 November 8 


F675W 


6714 


0.40 ±0.01 


1994 November 8 



Note. - Central wavelengths A c and observed 
(anomalous) flux ratios F f, g between images A2 and 
Ai in MG 0414+0534, in each of nine filters. The 2007 
November 3 observations were taken with the IMACS 
and PANI C instruments on the Magellan 6.5-m Baade 
telescope l IBate et al.ll2008T) . The 1997 August 14 ob- 
servations were taken with the NICMOS instrument on 
HST (obtained from the CASTLES Survey web page. 
The 1994 November 8 observati ons were taken with the 
WFP C2 instrument on HST JFalco, Lehar. & Shapiro] 
[T997D . 

scope. These data were first presented in IBate et all 
((2001) . 

2.2. SDSS J 0924+0219 

SDSS J0924+0219 is the most anomalous lensed quasar 
currently known. The minimum image A has been ob- 
served to be a factor of ~ 20 brig hter than the sad - 
dle point image D in the optical (Kee ton et al.l [2 006). 
The quasar was discovered by llnada et alj (|2003l ) in 
Sloan Digital Sky Survey (SDSS) imaging, and con- 
sists of an early-type lensing galaxy at zi — 0.394 
(lEigenbrod et a l. 2006) and a background quasar at z s — 
1.524 (jlnada et al.ll2003D . 

Again, we use three epochs of observational data. 
These are presented in Table d The 2008 March 
21 data were obtained by us using the Magellan 6.5 - 
metre Baade telescope (jFlovd. Bate, k Websterl [2009). 
The 2003 November 18-23 data were taken using the 
HST/NICMOS and W FPC2 instruments as part of 
the CASTLES Survey (IKeeton et all f2006h . T he 2001 
December 15 were obtained by llnada" "et~aTl ((2003) 
using the MagIC instrument on the Baade tele- 
scope, and re-reduced by us (details can be found in 
IFlovd. Bate, k Websterll2009l) . 

2.3. Q2237+0305 

Q2237+0305 is perhaps the most well-st udied gravita- 
tionally lensed quasar. It was discovered bv lHuchra et al.1 
(1985), and consists of a lensing galaxy at Zi = 0.0394 
and a background quasar at z s = 1.695. The two previ- 
ous sources had early type lensing galaxies; the lens in 
Q2237+0305 is a barred spiral. Near-perfect alignment 
between observer, lens and quasar results in four virtu- 
ally symmetric images of the background source, located 
in the bulge of the lensing galaxy. The optical depth to 
stars is therefore quite high, making the system an ex- 
cellent target for microlensing analyses. Typically, this 
is taken to mean the smooth matter percentage in mi- 
crolensing simulations can be set to zero. We will test 
this assumption here. 
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Table 2 

Observed flux ratios in SDSS J0924+0219 



Filter 


A c (A) 


Fobs = D/A 


Date 


u 


1 fi^nn + 1 450 


n OS + 05 


2008 March 21 


j 


i o^nn 4- son 


n 1 5 4- o 05 


9008 March 91 


Y 


1 (1200 + 500 


n 14 4- 05 


2008 March 21 


z' 


9134 ± 800 


0.19 ±0.10 


2008 March 21 


%' 


7625 ± 650 


0.16 ±0.10 


2008 March 21 


r' 


6231 ± 650 


0.10 ±0.10 


2008 March 21 


!/ 


4750 ± 750 


0.08 ± 0.08 


2008 March 21 


i' 


7625 ± 650 


0.08 ±0.05 


2001 December 15 


r' 


6231 ± 650 


0.07 ±0.05 


2001 December 15 


g' 


4750 ± 750 


0.06 ±0.05 


2001 December 15 


u' 


3540 ± 310 


<0.09 


2001 December 15 


F160W 


15950 ± 2000 


0.08 ±0.01 


2003 November 18 


F814W 


8269 ± 850 


0.05 ±0.005 


2003 November 18-19 


F555W 


5202 ± 600 


0.05 ±0.01 


2003 November 23 



Note. — Central wavelengths A c and observed (anomalous) 
flux ratios F obs between images A and D in SDSS J0924+0219, 
in each filter. The 2008 March 21 observations were taken with 
the IMACS and PANIC instruments on the Magellan 6.5-m 
Baade telescope .Floyd, Bate. & Webster1[2009f) . T he 2001 De- 
cemb er 15 data were taken using MagIC on Baade lllnada et al.l 
2003). The 2003 data were take n using HST/NICM OS and 
WFPC2 for the CASTLES Survey (Keeton et al.ll2006l 1. 

Q2237+0305 also differs from the two previous sources 
in that it does not contain a close image pair displaying 
a flux ratio anomaly. Nevertheless, its images are known 
to vary in brightness due to microlensing. We choose to 
examine images A and B both because they are well mod- 
elled, and they are roughly equi-distant from the centre 
of the lensing galaxy. 

Our observa ti onal data were obtained from 
lEigenbrod et al.l (|2008t) . in which 43 epochs of 
Q2237+0305 spectroscopic data from the FORS1 
spectrograph on the Very Large Telescope (VLT) were 
analysed. Eigenbrod and collaborators deconvolved 
their Q2237+0305 spectra into a broad emission line 
component, a continuum emission component, and 
an iron pseudo-continuum. The continuum emission 
was fit with a power-law of the form /„ cx u a " . The 
power-law fit was then split into six wavelength bands, 
each with a width of 250A in the quasar rest frame, and 
integrated in each band. The result is pseudo-broadband 
photometry in six wavebands, with contamination from 
broad emission lines and the iron continuum removed. 
We selected two epochs from these data for our analysis, 
separated by approximately a year. T he flux ratios are 
presen ted in Table [3] Following the lEigenbrod et al.l 
(2008) numbering, the 2005 November 11 dataset is 
epoch number 17, and the 2006 November 10 dataset is 
epoch number 28. 

3. MICROLENSING SIMULATIONS 

The simulation technique used in thi s work has 
been presented previou s ly in iBate et ail (|2008f ) and 
IFlovd. Bate, fc Websterl (|2009f ). In those papers, we 
marginalised over the smooth matter percentage s as a 
nuisance parameter. Here, we instead marginalise over 
the accretion disc radius <jq and the power-law index £ 
relating the radius of the accretion disc to the observed 
wavelength. 

Our microlensing simulations were con- 
ducted using an inverse ray-shooting tech- 



nique (IKavser. Refsdal fc Stabelll 1986; 

iWambsganss. Paczvnski fc Katzl I1990D . The key 
parameters in these simulations are the convergence 
Ktat, which is a measure of the focussing power of the 
lens, and shear 7, which is a measure of the distortion 
introduced by the lens. The lensing parameters used in 
this analysis are presented in Table SJ 

The convergence can be split into two components 
K-tot — + k>si where k* describes a clumpy stellar com- 
ponent and smoothly distributed component. We 
define the smooth matter percentage s to be the ratio 
of the continuously distributed component to the total 
convergence: 



fl 



We allowed the smooth matter percentage to vary from 
to 99 per cent, in 10 per cent increments. The smooth 
matter percentage is thus relatively coarsely sampled; 
our simulations were optimised to probe the accretion 
disc sources in each system. 

The microlenses in our simulations were drawn from a 
Salpeter mass function dN/dM cx M~ 2 35 with a mass 
range M max / M m i n = 50. Physical sizes are therefore 
scaled by the average Einstein Radius projected on to 
the source plane 770, which varies from system to system. 
Magnification maps were generated covering an area of 
24770 x 24?7o, with a resolution of 2048 x 2048 pixels. Ten 
maps were generated for each image and smooth matter 
percentage. 

As discussed in IBate et ail (|2008h . we randomly se- 
lected source positions in each combination of magnifi- 
cation maps, to build up a simulated library of flux ratio 
curves as a function of wavelength. Comparing these 
with the observations allows us to construct a three di- 
mensional likelihood distribution for the observed flux 
ratio spectrum F obs given three model parameters: the 
radius of the quasar source in the bluest filter o~o, the 
power-law index relating observed wavelength to radius 
of the source and the smooth matter percentage s. We 
can convert these likelihoods to an a posteriori probabil- 
ity distribution for the three model parameters given the 
observations using Bayes theorem: 



d 3 P 



doodCds 



ex L(F° bs \aoX, s) dP / lio1 dP P; ior dP P rior (2) 



dor. 



d( 



ds 



Uniform priors were used for the two dimensionless 
quantities, smooth matter percentage s and power-law 
index A logarithmic prior was used for the radius 
ctq . W e note that this d iffers slightly from the anal- 
yses i n IBate et al.l (|2008D and IFlovd. Bate, fc Websterl 
(2009), where a uniform prior was also used for a . We 
will briefly discuss prior dependence in Section |H We 
marginalise over the accretion disc parameters o-q and £ 
to obtain a probability distribution for the smooth mat- 
ter percentage s: 



dP 
ds" 



d 3 P 



doodCds 



der dC 



(3) 



Our analysis focusses on two lensed images in each sys- 
tem. By dealing with flux ratios between images only, we 
remove the intrinsic quasar flux from the problem, pro- 
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Table 3 

Observed flux ratios in Q2237+0305 



Band Emitted A c (A) Observed A c (A) B/A Date 



1 1625 ± 125 4379 ± 337 0.52 ± 0.02 2005 November 11 

2 1875 ± 125 5053 ± 337 0.51 ± 0.02 2005 November 11 

3 2125 ± 125 5727 ± 337 0.50 ± 0.01 2005 November 11 

4 2375 ± 125 6401 ± 337 0.49 ± 0.01 2005 November 11 

5 2625 ± 125 7074 ± 337 0.48 ± 0.01 2005 November 11 

6 2875 ± 125 7748 ± 337 0.47 ± 0.01 2005 November 11 

1 1625 ± 125 4379 ± 337 0.33 ± 0.02 2006 November 10 

2 1875 ± 125 5053 ± 337 0.34 ± 0.02 2006 November 10 

3 2125 ± 125 5725 ± 337 0.35 ± 0.01 2006 November 10 

4 2375 ± 125 6401 ± 337 0.36 ± 0.02 2006 November 10 

5 2625 ± 125 7074 ± 337 0.37 ± 0.02 2006 November 10 

6 2875 ± 125 7748 ± 337 0.37 ± 0.01 2006 November 10 



Note. — Two epochs of observ ational B/A flux ratios for Q2237+0305, ob- 
tained from Eigcnbrod ct al. (2008). Photometry was extracted from spectra 
obtained with the FORS1 spectrograph on the Very Large Telescope (VLT) at 
the European Southern Observatory (ESO). Filter wavelengths are in the rest 
frame of the source qu asar, located at a redshift of z s = 1.695. Following the 
Eigcnbrod ct al. (2008) numbering, the 2005 November 11 datasct is epoch num- 
ber 17, and the 2006 November 10 dataset is epoch number 28. 



Table 4 

Lcnsing parameters for the images of interest in 
this analysis 



Quasar 


Image 


Ktot 


7 


fltot 


MG 0414+0534 


A! 


0.472 


0.488 


24.2 


MG 0414+0534 


A-2 


0.485 


0.550 


-26.8 


SDSS J0924+0219 


A 


0.502 


0.458 


26.2 


SDSS J0924+0219 


D 


0.476 


0.565 


-22.4 


Q2237+0305 


A 


0.413 


0.382 


5.03 


Q2237+0305 


B 


0.410 


0.384 


4.98 



Note. — Convergence Ktot, shear 7 and magnifi- 
cation fitot for each of the lensed images examined in 
this paper. A negative total magnification is inter- 
preted as a parity flip. MG 0414+0534 par a meter s 
were obtained from Witt. Mao & Schcchtcr (1995). 
SPSS J0924+0219 parameters were obtained from 
IKeeton et aIU2006l ). Q2237+0305 were obtained via 
private commu nication with Cath ryn Trott, based 
on modelling in lTrott etaLl lf2009T ). 

vided the difference in light travel time between images is 
short. We assume that the smooth matter percentages in 
the two lensed images are identical. This is reasonable for 
the two anomalous systems, MG 0414+0534 and SDSS 
J0924+0219, as the anomalous images lie very close to 
each other. In Q2237+0305, where the images are widely 
separated, we choose to analyse image A and B only as 
they are essentially equi-distant from the centre of the 
lensing galaxy, and do not lie atop any obvious spiral 
features. 

The probability distributions we obtain for smooth 
matter percentage are presented for MG 0414+0534 (Fig- 
ure [J), SDSS J0924+0219 (Figure [2]) and Q2237+0305 
(Figure [3]). The dashed histograms show the differential 
probability distributions, and the solid lines show the 
cumulative probability distributions. 

4. RESULTS AND DISCUSSION 

We obtain the following formal constraints on smooth 
matter percentage at the image positions in each system: 




+10 



501^ per cent in MG 0414+0534, 80 
J0924+0219, and < 50 in Q2237+0305 (68 per cent con- 



10 per cent in SDSS 



smooth matter percentage 

Figure 1. Probability distribution for smooth matter percentage 
in MG 0414+0534. The differential distribution (dashed line) and 
cumulative distribution (solid line) are both provided. 



fidence limits are quoted). Our simulations are not op- 
timised to probe smooth matter percentage; we sample 
smooth matter parameter space only sparsely, in order to 
reduce simulation time. These results should therefore be 
considered estimates, rather than exact measurements. 
Nevertheless, they provide an interesting, and currently 
poorly explored, measurement of smooth matter content 
within only a few effective radii of early type galaxies. 

In MG 0414+0534, the differential probability distribu- 
tion (Figure [TJ dashed line) does not particularly favour 
any single smooth matter percentage. This leads to a 
quite broad formal constraint on the smooth matter per- 
centage. 

Conversely, our measured smooth matter percentage in 
SDSS J0924+0219 is high, as we would expect for such 
an anomalous system. There is one other measurement 
of sm ooth matter percent age of this system in the litera- 
ture. IKeeton et al.l (|2006l ) estimated it to be 80 to 85 per 
cent, based on the estimated size of the broad emission 
line region in the system. Our analysis is focussed on the 
accretion disc, and finds essent ially identical results. We 
do use the IKeeton et al.l (|2006l ) observational data in ob- 
taining our constraints, however excluding it and working 
only with the Magellan data does not significantly alter 
our results. 
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smooth matter percentage 

Figure 2. Probability distribution for smooth matter percentage 
in SDSS J0924+0219. The differential distribution (dashed line) 
and cumulative distribution (solid line) are both provided. 



age positions. To do this, we take the observed effective 
radius of the lensing galaxy R e , a nd compare it with 
Figure 10 in iBernardi et al.l (|2003t ) to obtain a g-band 
surface brightness in magnitudes per square arcsecond 
(assuming no evolution in early type galaxies between 
r edshift z = and t h e reds hift of the lensing galaxy) . 

Kauffma nn et al.l (|2003|) provide a relationship be- 
tween mass-to- light ratio M/L and <7-band magnitude 
derived from 10 5 Sloan galaxies in their Figure 14. 
Using the q-band s urface brightness obtained from 
IBernardi et al.l (|2003D . we can get a rough mass-to-light 
ratio for our lensing galaxy. We convert our surface 
brightness from magnitudes per square arcsecond into 
solar luminosities per square parsec using the following 
relationship: 




smooth matter percentage 

Figure 3. Probability distribution for smooth matter percentage 
in Q2237+0305 at the image A and B positions. The B/A flux ratio 
was used as those images are at roughly the same projected distance 
from the centre of the lensing galaxy (~ 0'.'95). The differential 
distribution (dashed line) and cumulative distribution (solid line) 
are both provided. 



As has been discussed earlier, the lensed images in 
Q2237+0305 lie in the bulge of the lensing galaxy. Stars 
are therefore expected to dominate the microlensing sig- 
nal, rather than a smooth matter component. Indeed, 
we find a smooth matter percentage that is consistent 
with zero in this system. There is a peak in the differ- 
ential probability distribution (Figure 02 dashed line) at 
~ 20 per cent. Though we are reluctant to suggest that 
this peak is real, we do note that such a feature could be 
evidence of additional a bsorbing features along the line 
of sight (see for example iFoltz et al.lll992l who reported 
Mgll absorption features in the spectrum of Q2237+0305 
at a redshift of 0.97). 

To confirm that our results are not dominated by the 
choice of prior probability for the radius of the back- 
ground quasar accretion disc a$, we repeated our analy- 
sis using a uniform prior rather than a logarithmic prior. 
Within their errors, we found no significant variation in 
our constraints on smooth matter percentage in any of 
our systems. 

We can perform a simple calculation to obtain a rough 
theoretical prediction of the percentage of dark matter we 
expect to see at the image pos itions in each source. The 
method is briefly described by Schechter & Wambsganss 
(2002), but will be repeated here for clarity. It is only 
applicable to MG 0414+0534 and SDSS J0924+0219 as 
it assumes an early type lensing galaxy. We begin by 
working out a stellar surface mass density S s at the im- 



^[mag/arcsec 2 ] = M Qg + 21. 572-2. 51og lo S[L /pc 2 ] (4) 

where S is the surface brightness and M Qs = 5.45 is the 
solar magnitude in the g-band. Now that we have the 
lens galaxy surface brightness in solar magnitudes, we 
can use the mass-to-light ratio to convert into a stellar 
surface mass density in solar masses per square parsec. 
This is the stellar surface mass density at the effective 
radius R e - we use the deVaucoulers profile to extrapo- 
late this value out to the Einstein Radius, which is the 
approximate location of the lensed images. 

The final piece of information we need is the surface 
mass density at the image positions. The critical sur- 
face mass density E cr for lensing is obtained from the 
following relationship: 



4ttG D d D ds 



(5) 



where D d is the angular diameter distance to the lens, 
D s is the angular diameter distance to the source, and 
Dd s is the angular diameter distance between lens and 
source. The other symbols have their usual meanings. 
For an isothermal sphere, the surface mass density at the 
Einstein Radius is half the critical surface mass density 
E cr . The smooth matter percentage s is simply obtained 
as follows: 



1 



0.5Scr 



(6) 



MG 041 4+0534 has a deVauc oulers effective radius 
R e = 0"78 (|Kochanek et al.|[2000h . which gives a g-band 
magnitude of ~ 21.75 and a mass-to- light ratio Mj L ~ 4. 
This gives a stellar surface mass density at the effective 
radius of 514M^/joc 2 . For an Einstein Radius of l'/15 
(jTrotter. Winn fcH cwitt 200(3), the stellar surface mass 
density on the Einstein Ring is smaller by a factor of 0.46, 
giving E s = 236M Q /pc 2 . For a source at z = 2.64 and 
a lens at z = 0.96, the critical surface mass density for 
lensing is S cr = 2189Af Q /pc 2 . This gives a theoretically 
predicted smooth matter percentage of 78 per cent. Note 
th at this number differs slightly fro m the figure quoted 
in iSchechter fc Wambsganss) (120021) a s they were using 
prepri nts of IBernardi et alT~(|2003[ ) and Kauffmann et alj 
( 2003). This prediction is on the high side of our mea- 
sured smooth matter percentage of 50^4° per cent, al- 
though it is formally consistent. 
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SDSS J0924+0219 has two meas ured effective radii in 
the literature: R e = 0"3 1 ± 0"02 ([Morgan et al.l l200l 
and R e = 0"50 ± 0'.'05 (lEigenbrod et al.1 120061). Both 
result s were obtained using HST data. iMorgan et al.1 
(2006) fit a deVaucouler s profi le to the lensing galaxy, 
whereas lEigenbrod et al.l ((2006) chose a two-dimensional 
exponential disk. Since an exponential profile is shal- 
lower than a deVaucoulers profile, it would tend to give 
a larger effective radius. This may allow the two effective 
radii to be reconciled, however we will deal with each case 
separately. The system has a lensing galaxy at z = 0.39 
and a source at z = 1.524, giving a critical surface mass 
density of S cr = 2323Af /pc 2 . 

For an effective radius R e = 0'.'50 we expect a g-band 
surface brightness of ~ 20.5mag/arcsec 2 , and a mass-to- 
light ratio M/L — 5.6. At the effective radius, the stellar 
surface mass density is 2285M P ) /pc 2 . For an Einstein Ra- 
dius of 0'.'85 (jlnada et al.ll2003[ ) and a deVaucoulers pro- 
file, the stellar surface mass density at the Einstein Ra- 
dius is lower by a factor of 0.35, giving S s = 791M Q /pc 2 . 
Thus, we predict a smooth matter percentage of 32 per 
cent. This is significantly lower than our la measured 
value of 80^10- Even at the 95 per cent level, our mea- 
sured smooth matter percentage is not consistent with 
the prediction of this rough calculation. 

The situation is somewhat different for an effective ra- 
dius R e = 0'.'31. Here, an early type galaxy with no 
evolution between z — and z = 0.39 should have a 
g-band surface brightness of ~ 19.75mag/arcsec 2 and a 
mass-to-light ratio M/L — 5. This gives a stellar sur- 
face mass density of 4052M Q /pc 2 at the effective radius. 
Again assuming a deVaucoulers profile, the stellar sur- 
face mass density is reduced by a factor of 0.13 at the 
Einstein Radius of 0'.'85, giving 539M©/pc 2 . We there- 
fore expect 46 per cent of the mass to be in stars, and 54 
per cent in smooth matter. This is consistent with our 
measured smooth matter percentage at the 95 per cent 
level (80±£°). 

In Figure |4] we show a comparison between the results 
of this rough theoretical calculation and our microlens- 
ing measurements for MG 0414+0534 (circle) and SDSS 
J0924+0219 (square). No errors are provided for the 
rough theoretical calculations, although we emphasise 
that they should only be considered estimates. With only 
two data points, there does not seem to be any system- 
atic difference between the microlensing measurements 
and the rough calculation. 

It is, however, not correct to connect smooth matter 
percentage in our simulations directly with dark matter 
content in the lens. Intervening systems (which may not 
be detected otherwise) can contribute smooth matter to 
the surface mass density, as can gas and dust in the lens. 
As the lensing galaxies in both MG 0414+0534 and SDSS 
J0924+0219 are early type galaxies, we would expect 
this contribution to be ve ry small. More importantly, 
iLewis &; Gil-Merino (2006) showed that small compact 
masses can mimic a smooth matt er component. This wa s 
first demonstrated analytically in lDyer fc Roederl (119811). 
For a source with a radius 0.1%) ILewis fc Gil-Merinol 
(2006) find that compact masses smaller than ~ 0.01M Q 
can mimic a smoothly distributed mass component. This 
may help to explain why our measured microlensing 
smooth matter percentage in SDSS J0924+0219 appears 




20 40 60 80 100 
microlensing 

Figure 4. Comparison between smooth matter percentages ob- 
tained from microlensing measurements and rough theoretical cal- 
culations in this paper. No errors are given for the theoretical 
calculations (see Section [4j, although they should only be consid- 
ered rough estimates. The dashed line shows 1:1 correspondence. 
Results for MG 0414+0534 (circle) and SDSS J0924+0219 (square) 
are shown, with 68 per cent confidence limits. 

to be slightly higher than we would expect from a simple 
deVaucoulers galaxy profile in an isothermal dark matter 
profile. 

We note that the three lensing systems analysed in this 
paper do not represent a fair statistical sample of gravita- 
tionally lensed quasars. The single-epoch imaging tech- 
nique used here requires short time delays between im- 
ages, and so we preferentially select lenses with close im- 
age pairs (or high symmetry in the case of Q2237+0305). 
We have chosen close image pairs that exhibit strongly 
anomalous flux ratios for our analysis. Since a priori, 
anomalo us flux ratios are driven by a smooth matter 
fraction (|Schechter fc W ambsganss 2002), this selection 
biases our estimates towards larger smooth dark matter 
fractions. Thus our results cannot be used to make global 
statements regarding smooth dark matter in galaxies. 
However eight out of ten close image s ystems display 
anomalous flux ratios (|Poolev et al.l|200 7). and so we ex- 
pect the smooth matter percentage at those image po- 
sitions to be high. This indicates that any bias in our 
analysis of individual systems would be relatively minor. 
In the future our analysis could be extended to a statis- 
tical sample through the use of m onitoring data an d the 
light curve technique discussed in iKocha nck (2003). 

5. CONCLUSIONS 

We have presented estimates of the smooth matter per- 
centage in three lensing galaxies along the line of sight to 
the lensed images. We find a smooth matter percentage 
of 50±^ in MG 0414+0534, 80±^ in SDSS J0924+0219, 
and < 50 in Q2237+0305, with 68 per cent confidence. 
In the two systems where the lensed images lie in the 
outer regions of the lensing galaxies (5 to 10 kpc from 
their centres), these measurements are inconsistent with 
zero smooth dark matter in the lensing galaxies. In 
Q2237+0305, where the lensed images lie in the central 
bulge of the lensing galaxy and so stars are expected to 
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dominate the microlensing signal, our result is consistent 
with zero per cent smooth matter, as expected. 

These measurements were obtained using a single- 
epoch imaging technique that is free from the need for 
long-term monitoring campaigns, which are required for 
detailed analysis of the lens mass profile. Our results 
also do not depend upon the typically unknown veloci- 
ties of the stars in the lensing galaxy, and of the lensing 
galaxy and background source themselves, which enter 
into analyses of microlensing lightcurves. It is however 
only appropriate for systems with time delays between 
images of less than a day, so we can ensure we are ob- 
serving the background source in the same state along 
each line of sight. It also does not provide us with any 
information on the slope of the dark matter profile; for 
that, stellar velocity dispersions in the lensing galaxy are 
required. 

Nevertheless, our technique does provide an observa- 
tionally inexpensive method for estimating the dark mat- 
ter fraction in lensing galaxies at the location of lensed 
images. Gravitational lensing remains the only method 
for directly probing the dark matter content of galaxies. 
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